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ABSTRACT 

We investigate how cluster morphology is affected by the cosmological con- 
stant in low-density universes. Using high- resolution cosmological N-body/SPH 
simulations of flat (fio = 0.3, Aq = 0.7, ACDM) and open {Qq = 0.3, Aq = 0, 
OCDM) cold dark matter universes, we calculate statistical indicators to quan- 
tify the irregularity of the cluster morphologies. We study axial ratios, center 
shifts, cluster dumpiness, and multipole moment power ratios as indicators for 
the simulated clusters at z = and 0.5. Some of these indicators are calculated 
for both the X-ray surface brightness and projected mass distributions. In ACDM 
all these indicators tend to be larger than those in OCDM at z = 0. This result 
is consistent with the analytical prediction of Richstone, Loeb, & Turner, that is, 
clusters in ACDM are formed later than in OCDM, and have more substructure 
at z = 0. We make a Kolmogorov-Smirnov test on each indicator for these two 
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models. We then find that the results for the multipole moment power ratios 
and the center shifts for the X-ray surface brightness are under the significance 
level (5%). We results also show that these two cosmological models can be 
distinguished more clearly at 2; = than z = 0.5 by these indicators. 

Subject headings: galaxies:clusters:general - cosmology:theory - methods:numerical 

1. INTRODUCTION 

Measurement of the cosmological parameters, such as the Hubble constant Hq, the den- 
sity parameter Qq, and the cosmological constant A, is one of the most important studies 
in observational cosmology. From observational evidence, the density parameter Qq is es- 
timated to be around 0.3 (e.g. Bahcall et al. 1999). The inflationary cosmology requires 
^^0 + -^0 = 1; where Aq = A/ (Sifo), and is supported by recent observations of the cosmic mi- 
crowave background (e.g. de Bernardis et al. 2000). Distant Type la supernovae (SNe la) are 
assumed to be standard candles to explore the effect of a cosmological constant by searching 
for evidence of accelerated expansion (Perlmutter et al. 1999; Riess et al. 1998). Several 
systematic uncertainties in the observational data of SNe la have been pointed out (Totani 
& Kobayashi 1999; Gibson et al. 2000). Therefore, independent studies in the non-linear 
regime are important to confirm the values of the cosmological parameters. 

Richstone, Loeb, & Turner (1992) proposed that the fraction of clusters which have 
significant substructure can be a probe of the density parameter, Qq, since the fraction 
of recently formed clusters strongly depends upon Oq ^-nd such young clusters are likely 
to have substructure. They have also revealed that the formation epoch of clusters in a 
universe with a cosmological constant is later than that in a universe without a cosmological 
constant, although the effect of the cosmological constant is much smaller than that of the 
density parameter. In fact, Wilson, Cole, & Prenk (1996) showed that the statistics of the 
quadruples of the column mass density of clusters strongly depends on the value of JIq, while 
it is quite insensitive to the value of Aq. 

Some authors have studied if cosmological models with and without a cosmological con- 
stant can be discriminated by the statistics of the irregularity of cluster morphologies (Jing 
et al. 1995; Crone, Evrard, & Richstone 1996; Buote & Xu 1997). They performed pure 
N-body simulations in low-density ffat universes (Qq — 0.2-0.35, Aq = 1 — ^0) and open 
universes {Qq = 0.2-0.35, Aq = 0). Jing et al. (1995) and Crone et al. (1996) quantified 
morphologies of X-ray surface brightness of clusters from the A^-body results assuming hy- 
drostatic equilibrium of the hot gas in the gravitational potentials of the clusters. Using 
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a Kolmogorov-Smirnov test (hereafter KS-test), they showed that a low-density flat uni- 
verse and an open universe are distinguishable from each other by the center shifts. On 
the other hand, Buote & Xu (1997) obtained multipole moment power ratios for clusters 
in their N-body simulations and compared those in a low-density flat universe and an open 
universe. They concluded that these two cosmological models are not distinguishable by the 
multipole power ratios. Since these authors estimated the hot gas distribution from N-body 
simulations, their results are in question, and then simulations including hydrodynamics are 
desired. 

In previous hydrodynamic simulations (Evrard et al. 1993; Mohr ct al. 1995) the num- 
bers of clusters were not enough to perform a statistical comparison. Although Valdarnini, 
Ghizzardi, & Bonometto (1999) analyzed a sufficient number of clusters for various cosmo- 
logical models, the effect of the tidal field from larger scales (> 30-50/i'"^Mpc) was neglected 
in their simulations because they performed tree-SPH simulations of individual clusters se- 
lected from the large A^-body simulations. By comparing the simulated chisters and the 
ROSAT X-ray cluster data, they pointed out that the statistical test for the observed data 
favors the high-density (Qq = 1, SCDM) model rather than the concordant low-density flat 
{Qq = 0.3, Ao = 0.7) model. It is worthwhile to conflrm their results by simulations that take 
account of the tidal fleld from outside the clusters. We then simulate comoving 150 h^^Mpc 
boxes with particle-particle-particle-mesh (P'^M)-Smoothed Particle Hydrodynamics (SPH) 
code. 

In this paper, we study morphologies of simulated clusters in flat {Qq = 0.3, Aq = 
0.7, ACDM) and open {flo = 0.3, Aq = 0, OCDM) cold dark matter universes by using 
high resolution simulations. We perform large simulations to obtain a sufficient number 
of clusters for statistical comparisons. A reason why we study the difference between low- 
density universes despite the fact that recent CMB data (e.g. de Bernardis et al. 2000) 
supports a flat universe is as follows. As we mentioned above, the irregularity of cluster 
morphologies is quite sensitive to the value of Qq while it is rather insensitive to the value 
of Aq. Then the statistical indicators that can discriminate between the low-density models 
with and without the cosmological constant Aq enable us to detect the small difference 
in the formation histories of clusters in different but similar cosmologies such as the low- 
density universe with a time dependent vacuum energy. Furthermore, because the CMB data 
provide us only the information about the linear density fluctuations, the consistency test 
in the non-linear regime is very important. For instance, if the flnding of Valdarnini et al. 
(1999) is true, we should consider what causes the discrepancy. Unfortunately, the resolution 
in our simulations is not sufficient to make a direct comparison with high-resolution X- 
ray observations by Chandra or XMM and our simulations do not include some important 
physical processes like the cooling and heating, which may affect the morphology of the 
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X-ray surface brightness. Therefore here we restrict ourselves to finding effective statistical 
indicators to discriminate between two low-density universes and we do not compare our 
simulation results with the observations. 

We identify clusters in our simulation and calculate various statistical indicators for 
them. The KS-tests are performed to measure how effectively these indicators distinguish 
between these two cosmological models. These indicators are calculated for the projected 
mass density as well as the X-ray surface brightness of clusters of galaxies. The X-ray 
surface brightness reflects the distribution of the hot gas in the cluster. On the other hand, 
the projected mass density is dominated by the distribution of dark matter. Although dark 
matter cannot be observed directly, a method by which the mass distribution is reconstructed 
from small distortion images of background galaxies caused by gravitational lensing (these 
small distortions are called the 'weak shear field') has been developed recently (e.g. MeUier 
1999; Luppino & Kaiser 1997; Clowe et al. 2000). The precision of this method is not yet 
sufficient to probe cluster substructure. It is, however, likely that this mass reconstruction 
method will progress and become a useful tool to measure substructure in galaxy clusters. 
Schneider & Bartelmann (1997) have already investigated the possibility of deriving multipole 
moments of the projected mass distribution of clusters with weak-lensing. 

The organization of this paper is as follows. In §2, we describe the numerical simulations 
used in this paper, the method of cluster identification, and the definitions of indicators which 
we use. We give mean values and standard deviations of the indicators for simulated chistcrs 
and show results of the KS-tests on these indicators in §3. In §4, we discuss our results and 
present our conclusion. 

2. METHOD 

2.1. Numerical Simulations 

Our simulations are based on a P^M-SPH algorithm. Detailed description of our simu- 
lation method is given in Yoshikawa, Jing, & Suto (2000). 

We use two cosmological models, OCDM and ACDM. Cosmological parameters in the 
two models, except for the cosmological constant, are as follows: the Hubble constant in 
units of 100km s~^ Mpc^^, h — 0.7, the density parameter, Qq = 0.3, the baryon density 
parameter, Qf, — 0.015/i""^, the rms density fluctuation amplitude on a scale 8/i~^Mpc, 
CTg = 1.0, and the power-law index of the primordial density fluctuations, n — 1.0. The 
normalized cosmological constant, Aq is and 0.7 for OCDM and ACDM, respectively. The 
simulation of ACDM is the same as L150A in Yoshikawa et al. (2000) and that of OCDM is 
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carried out for this paper. 

Each simulation employs A^dm = 128^ dark matter particles and the same number of 
SPH particles. The mass of a dark matter particle and an SPH particle are 1.7 x 10^^ Mq and 
2.0 X IO^^Mq, respectively. The size of comoving simulation box, Lbox, is 150/i~^Mpc, and 
the box has the periodic boundary conditions. We use the sphne (S2) gravitational softening 

(Hockney & Eastwood 1981), and the softening length, Cgi-av, is set to be Lbox/(10-/V^^) 
(~ 120/i~^kpc). The smoothing length of each SPH particle is determined by 

where hf^^ is the smoothing length of the i-th particle at the n-th time step, A^^*^"''' is the 
number of neighbor particles of the i-th particle at the n-th time step (the number of particles 
inside a sphere of radius h[^''), and A^^ is mean number of neighbor particles (Hernquist & 
Katz 1989). We set Ng = 32 and the minimum SPH smoothing length as /imin = egrav/4 (~ 
30/i~^kpc). We use the COSMICS package (Bertschinger 1995) to generate initial conditions 
at 2; = 25. Our simulations are carried out on a VPP5000 and each run took about 50 hours 
using 8PE with parallel P^M-SPH. The memory resource needed to run each simulation is 
500MB. 

We calculate statistical indicators for the clusters in our simulations at 2; = 0.5 and 0. 
The reason for this is that the formation rate of galaxy clusters in ACDM is expected to 
exceed that in OCDM for z < 0.8 as shown in Fig.l which shows the formation rate of galaxy 
clusters as a function of redshift according to the analytical formula given by Richstone et 
al. (1992). In Fig.l the number of formed clusters is normalized by the present value and the 
cosmic time is also normalized by the present age of the universe. The solid and dashed curves 
indicate the formation rates in ACDM and in OCDM, respectively. Since just before z = 0.5 
the formation rate in ACDM exceeds that in OCDM, we expect that the difference due to the 
cosmological constant can already be detected at 2; = 0.5. Moreover, mass reconstruction 
based on observations of the week shear field has an advantage for high-redshift clusters 
rather than nearby clusters in principle. 

2.2. Cluster Identification 

The way of identifying galaxy clusters is as follows (Jing & Fang 1994; Thomas et al. 
1998). 

1. We pick up SPH particles whose gas densities are more than 200 times the background 
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1.5 1 
redshift 







Fig. 1. — The formation rate of galaxy clusters. The sohd curve indicates the formation rate 
in ACDM and the dashed curve indicates that in OCDM. 
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(gas and dark matter) density. Here, we define 'gas density of an SPH particle' as gas 
density at the position of the SPH particle. Processes at step 2 and 3 are performed 
for these selected particles. 

2. We perform a Priends-of-Priends method (hereafter PoP; e.g. Davis et al. 1985) for the 
dense SPH particles selected by step 1. The finking length, I, is defined as hn~^/^, where 
n is mean number density of all SPH and dark matter particles in the simulation box, 
and the constant parameter h = 0.5. It should be noted that the dense SPH particles 
are confined within small scale areas. Since the typical distance between these areas 
is significantly larger than / (= 6n^^/"^), the number of groups found by FoF does 
not depend on b strongly. We confirm that both the number of the groups found in 
this step and the number of clusters found in following steps are not so different for 
0.2 < 6 < 1.0. 

3. The densest SPH particle in each group is defined as a 'core' particle of the group. 

4. We draw a sphere of which the center is on the position of the 'core', and seek a radius 
in which the mean density of total matter is 200 times the background density of total 
matter. This radius is called r2oo- 

5. If the total mass in the sphere with radius r2oo is more than 2 x 10^^h~^MQ, a set of 
particles in the sphere is a candidate for being a cluster of galaxies. 

6. If another 'core' particle exists in the sphere, the set which belongs to the less dense 
'core' is removed from the cluster candidate list. 

The candidates which are not removed finally are identified as clusters of galaxies. We define 
the center of a cluster by the position of the 'core' particle. We find that this position is 
very close to the particle having the minimum gravitational potential in the cluster and the 
distance between the 'core' particle and the potential minimum is at most several ten kpc 
in our simulation. We also find that the position of the 'core' is close to the center of mass 
of the cluster except for some clusters which have large substructures. 

Pig.2 shows the mass distribution of the clusters. The upper panel describes the dis- 
tribution in ACDM and the lower one describes the distribution in OCDM. The solid line 
indicates the number of clusters at 2; = and the dashed line indicates the number at z = 0.5. 
The maximum masses of clusters are 2.7 x lO^^h'^MQ (ACDM, z = 0), 9.3 x lO^^h'^MQ 
(ACDM, z = 0.5), 2.7x lO^^/i-^M© (OCDM, z = 0), and 1.2 x lO^^h'^MQ (OCDM, z = 0.5). 
Table 1 shows the number of clusters found in each model. The numbers of clusters increase 
from z = 0.5 to z = 0. Table 2 is the ratio of the mass inside clusters to the total mass in 
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the simulation box. These ratios in both models are similar to each other at 2; = 0, but the 
ratio in ACDM is smaller than that in OCDM at z — 0.5. This is consistent with the fact 
that ACDM clusters form later than OCDM clusters. 



2.3. Calculation of Column Density and X-ray Surface Brightness 

Two of the indicators we use in this paper, the center shifts and the power ratios, are 
calculated for X-ray surface brightness and the projected mass density for simulated clusters. 
Our numerical projection method from 3D to 2D is as follows. 

At first, for each cluster, we prepare a cube centered on a core particle and the sides of 
the cube are set to 2r2oo- The total number of grid points in the cube is 128^. Each grid 
point is assigned indices k) and this point is expressed by rijk- We calculate gas density 
and temperature at each grid point from SPH particles in this cube. The gas density, Pgas, 
and temperature, T, at a grid point r^jk are given by 

PgUfijk) =^W{rs-ijk/hs)^, (2) 

T{r,,,) = V t^^}l^w{rs;jk/hs)us, (3) 
Kb 

s 

where hs, m^, Ug, and rs-^ijk are the s-th particle's smoothing length, mass, specific energy, 
and distance to the grid point {rs;ijk — \fijk ~ '''s], where is the position of the s-th 
particle), respectively. W{t) is an SPH kernel, which is defined as 

(3/2)t2 + (3/4)^3 if < t < 1 
tf/A if 1 < t < 2 

otherwise. 

We assume the specific heat ratio 7 = 5/3 and the mean molecular mass // = 0.6. The 
smoothing length of an SPH particle is calculated in the usual way of the SPH method 
(Hernquist & Katz 1989; Monaghan 1992). 



Model 


z=0 


z=0.5 


ACDM 


66 


28 


OCDM 


74 


38 




Table 1: The number of the clusters for each model and redshift. 
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ACDM 



- ACDM(z = 0) 
ACDM(z = 0.5) 



5x10'^ 10 
Moss of cluster(MQ/h) 

OCDM 



2x10 



OCDM(z = 0) 

OCDM(z = 0.5) 



5x10'* 10 
Moss of cluster(Mg/h) 



2x10 



Fig. 2. — Mass distribution of the clusters. The upper panel describes the distribution in 
ACDM and the lower one describes the distribution in OCDM. The solid line indicates the 
number of clusters at ^ = and the dashed line indicates the number at ^ = 0.5. 
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in order to obtain the mass density at each grid point from dark matter particles in each 
cluster, we use an interpolation technique as in the previous SPH method. The smoothing 
length of each dark matter particle is set in order that each dark matter particle has 32 ± 3 
neighbors. 

The column density, cr, is derived from the projection of the total mass density, Pdm+Pgas, 
where pdm is dark matter density and Pgas is gas density. The X-ray surface brightness. Ex, 
is derived from the projection of p^^T^^"^. 

2.4. Statistical Indicators 

We calculate statistical indicators for each cluster of galaxies to quantify substructure 
and perform the KS-tests. As statistical indicators, we adopt the axial ratio, the Mint, the 
multipole moment power ratio, and the center shift. Definitions are described below. 

2.4-1- Axial Ratio 

The axial ratio is an indicator to show deviation from sphericity of a cluster of galaxies 
(Dutta 1995; Jing et al. 1995; Thomas et al. 1998). In a coordinate system whose origin is 
at the center of mass of the cluster, the following tensor's eigenvalues are calculated for all 
dark matter particles and SPH particles. 

lij ^ ^ TflgXiXj . (4) 
s 

These eigenvalues are labeled Ai, A2, A3 in decreasing order. If a cluster is an ellipsoid, its 
axial ratio is obtained by (Ai/Aa)^/^. 



Model z=0 z=0.5 
ACDM 10.2% 3.64% 

OCDM 11.1'X 5.21'7o 



Table 2: The ratios of the mass involved in the clusters to total mass in each simulation. 
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24-2. Mint 

We use Mint as an indicator of dumpiness of clusters. We calculate Mint in a similar way 
to Thomas et al. (1998). For each cluster we perform an FoF method for all dark matter 
particles and SPH particles with initial linking length, / = Uc where Uc is the mean 
number density of dark matter and SPH particles in the cluster. Next I is gradually lowered 
until I becomes (lOOric)"^/^ This causes the cluster to break up into several subclumps. At 
each stage Mint is defined as 

Mint = , (oj 

mi 

where mi > m2 > m^ are masses of the three largest clumps. We use the maximum value 
of Mint 3S a nieasure of dumpiness of the cluster. 



2.4.3. Center Shift 

The center shift measures major deviations from symmetry in the cluster mass distri- 
bution (Jing et al. 1995; Crone et al. 1996). We calculate the center shift for each cluster 
in a slightly different way from Jing et al. (1995). We first search for a peak value, Cpeak, of 
a or Ex in the cluster. The lowest contour level, ciowest, is defined as the mean value of a 
or Ex at 0.5r2oo- Next the i-ih contour level is given as q = Cpeak(ciowest/cpeak)'''"'''°°*, where 
JT-cont is the total number of contours. 

The center shift, C, is defined as 

''^cont 

C= J] ^.{(x.-a;)^ + (!/,- l/)2}, (6) 
1=1 

where {xi,yi) is the center of the i-th contour, and x — Ylii'^i^ii V — T^i'^iUi- The 
weight of each contour, Wj, is proportional to the surface integral of a or Ex in the region 
between this contour and the adjacent outer contour. The center shift shows the emission- 
weighted dispersion of the centers of contours. If a cluster has a a lot of substructure, then 
the outer contour's center is expected to shift from {x, y) and the center shift becomes large. 



2.4.-4-- Power Ratio 

The power ratios quantify the shape of projected cluster potentials and are derived 
from their multipole expansions (e.g. Buote & Tsai 1995, 1996; Buote 1998; Valdarnini et 
al. 1999). 
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The two-dimensional potential, ^(i?, 0), and column density, a{R,(l)), are related by 
Poisson's equation, 

V2*(i?,0)=(7(i?,0), (7) 

where R and are projected polar coordinates about the center of mass of the cluster, and 
we ignore the constant factor in equation(7). The multipole expansion of \&(-R, 0) which 
relies on the distribution of the material interior to R (e.g. Buote & Tsai 1995) is 



/ 1 \ °° 1 
*(fl,^) = -a„l„(-]-i: — 

^ ' m=l 



mR"" 

X {am cos UKf) + bm sin m(f)), (8) 



where and are defined as follows: 



and 



[ a{x'){RT COSUKp' (fx' (9) 
Jr'<r 

bm= [ a{x'){R')"' sin m<j)' fx', (10) 
Jr'kr 

where x' = {R', </)'). The m-th moment power, P„i, is defined by integration of the square 
of the m-th term of the multipole expansion, equation (8), over the boundary of a circular 
aperture of radius Rap, 



for m > 0, and 



for 171 = 0. 



The moment power depends on not only the irregularity of the potential shape but also 
the magnitude of a. Thus we use the moment power ratio, Pm/Po, as an indicator of a 
cluster's irregularity. We define the unit of R, the radius in polar coordinates, as r2oo, so 
that the power ratio is independent of the size of the cluster. Since the origin of the polar 
coordinates is the center of mass of the cluster, the dipole moment. Pi, is vanished. The 
higher order (m > 4, in this paper) terms are affected by minor irregularities in cluster 
shapes, hence we calculate the power ratio for only m = 2, 3, and 4. Because the power 
ratios depend on Rap (Buote & Xu 1997; Valdarnini et al. 1999), we calculate the power 
ratios for some different values of Rap- 

We also use the same definition of the power ratio for Ex by replacing a with Ex in 
equations (9) and (10). 
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3. RESULTS 

According to the analytical results of Richstone et al. (1992), the typical formation epoch 
of galaxy clusters in ACDM is delayed to lower redshift than in OCDM. This delay clearly 
appears at low-z {z < 0.8-0.7, as shown in Fig.l). We then calculate the indicators described 
in §2 at 2; = and z — 0.5 and perform the KS-test for each set of indicators obtained for two 
cosmological models. The KS-test can be used as a statistical test to estimate the ability 
of an indicator to be used to discriminate between the model with and without Aq. The 
result of the KS-test is the probability of the null-hypothesis that two distributions of the 
indicator are generated from the same population (Press et al. 1988). The significance level 
of the KS-test adopted in this paper is 5%. If the result of the KS-test for an indicator is 
under this level, we regard that it is able to distinguish two cosmological models by using 
these indicator. 



3.1. Axial Ratio 

The mean values of the axial ratios are shown in Table 3. From z = 0.5 to 2; = the 
axial ratio decreases. At z = 0, the axial ratio in ACDM is larger than in OCDM. This 
suggests that the clusters in ACDM have more substructures than in OCDM as expected 
from the analytical prediction. 

The results of the KS-tests for axial ratios are 0.495 for z = and 0.661 for z = 0.5. 
Both are over the significance level(5%). Thus the axial ratio is not useful to distinguish 
these two cosmological models. 

Jing et al. (1995) performed KS-test for axial ratios, which were obtained in a slightly 
different way to us, and the result of their KS-test is under the significance level (0.13 x 10~^). 
In their calculations of principal axes and axial ratios, they take, at first, all particles within 
the virial radius of a cluster. They then calculate new principal axes and axial ratios using 
particles within an ellipsoid with the principal axes and axial ratios just determined. They 
repeat the same calculation for the updated ellipsoid until the axial ratios converge. We also 



redshift 


ACDM 


OCDM 


z^O 


1.608 


1.585 


z^O.5 


1.644 


1.696 



Table 3: The mean values of the axial ratio. 
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calculate axial ratios by their method and perform a KS-test on them. In contrast with their 
result, the result of our KS-test is over the significance level (0.150 for 2; = and 0.750 for 
z — 0.5) though our resolution is similar to theirs. The reason of this difference is unclear. 



The mean values of Mint are shown in Table 4. Mint in ACDM is larger than that in 
OCDM at z = 0. 

The results of the KS-tests for Mint are 0.770 for z = and 0.192 for z = 0.5. Both 
of them are over the significance level. Again we can not use this indicator to distinguish 
between these two cosmological models. 

The mean values and the results of KS-tests are consistent with those of Thomas et al. 
(1998). 



3.3. Center Shift 

The mean values and the standard deviations of the logarithm of the center shifts are 
shown in Table 5. The number of the contours, ncont, is varied from 4 to 8 to investigate the 
effect of the number of the contours. For large rZcont; i-e. a small interval of contour level, C 
is large. In the case of a small interval, the center shift easily shows small substructure near 
to the center of the cluster. 

The center shifts for cr, C^, are significantly larger than those for Ex, C'e^- The shapes 
of the Ex contours reflect the distribution of the gas. Gas in a cluster relaxes more quickly 
than the coUisionless dark matter (Prenk et al. 1999). This is the reason why the shapes 
of the Ex contours become rounder than those of a and the center shifts for Ex become 
smaller than those for a. 

The center shifts in ACDM are larger than those in OCDM at the same redshift. This 



redshift 


ACDM 


OCDM 


z = 


1.289 


1.188 


z = 0.5 


1.281 


1.189 



Table 4: The mean value of the Mint. 
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Model 


^cont 


log(a) 


log(C^s,) 


ACDM 








(^ = 0) 


4 


-2.65(±0.72) 


-3.43(±0.78) 




5 


-2.52(±0.69) 


-3.31(d=0.82) 




6 


-2.43(±0.67) 


-3.26(±0.80) 




7 


-2.38(±0.64) 


-3.22(±0.80) 




8 


-2.33(±0.63) 


-3.20(±0.77) 


OCDM 








(^ = 0) 


4 


-2.80(±0.70) 


-3.64(d=0.75) 




5 


-2.72(±0.66) 


-3.56(±0.77) 




6 


-2.62(±0.67) 


-3.49(±0.75) 




7 


-2.59(±0.63) 


-3.47(±0.76) 




8 


-2.54(d=0.62) 


-3.46(±0.76) 


ACDM 








{z = 0.5) 


4 


-2.87(±0.58) 


-3.73(±0.67) 




5 


-2.71(±0.60) 


-3.62(±0.64) 




6 


-2.67(±0.57) 


-3.53(±0.64) 




7 


-2.63(±0.55) 


-3.49(±0.67) 




8 


-2.59(±0.55) 


-3.45(±0.64) 


OCDM 








{z = 0.5) 


4 


-2.99(±0.68) 


-3.68(±0.77) 




5 


-2.89(±0.70) 


-3.59(±0.80) 




6 


-2.77(±0.68) 


-3.53(±0.76) 




7 


-2.73(±0.70) 


-3.46(±0.76) 




8 


-2.68(±0.69) 


-3.41(±0.80) 



Table 5: Mean values of logarithm of center shifts 



Note. — The first eolumn indicates cosmological model and redshift. The second cohimn shows the 
number of contours. The third and fourth columns are mean value of log(C) for the column density and the 
X-ray surface brightness, respectively. 
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cr(z = 0) * 

cr(z = 0.5) o 
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8 
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Fig. 3. — The results of the KS-tests for center shifts as a function of the number of contours. 
Asterisks indicate the resuhs at z = 0, and circles indicate the results at z — 0.5. The dotted 
line describes the significance level(5%) 
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result reflects that the formation epoch of clusters in ACDM is later than in OCDM. 

Fig.3 shows the result of the KS-tests for the center shifts. The upper panel describes 
the results for the center shifts for the column density and the lower one describes the results 
for the X-ray surface brightness. In Fig.3 the dotted hne indicates the signiflcance level(5%), 
asterisks indicate the results at 2; = 0, and circles indicate the results at z — 0.5. Using 
the center shifts for Ex at ^ = 0, we can distinguish the two cosmological models in this 
range of number of contours. The effect of contour level interval is not very signiflcant for 
the KS-test. 

In order to study the effect of the lowest level of contour, we also calculate the center 
shifts for Y^x and a in the two cases that the lowest contour level is the mean value at 0.4r2oo 
and 0.7r2oo- For the case of 0.7r2oo, the results of the KS-tests are over the significance level, 
except for one case with (z=0, ncont = 8)- For the case of 0.4r2oo, the results are similar 
to the case of the lowest level at 0.5r2oo- We conclude that the center shift within ~ 0.5r2oo 
is a useful tool to clarify the presence of a cosmological constant. 

3.4. Power Ratio 

We calculate the power ratios for various Rap (0.4, 0.5, 0.6, 0.8, and 1.0 times r2oo) in 
order to study how Rap affects the result of the KS-test. We show the mean values and the 
standard deviations of the logarithm of the power ratios in Tables 6 and 7, for the column 
density and the X-ray surface brightness, respectively. For small i?ap, the power ratios are 
large as shown in Table 6 and 7. The most likely explanation of this property is that the 
power ratios for small Rap readily reflect the small substructures around the cluster center. 

We also calculate the power ratios for constant i?ap, e.g. 0.4/i"^Mpc, and 0.6/i^^Mpc. 
As shown in Table 8, the power ratios for the constant Rap are smaller than the proportional 
Rap case because the constant Rap is less sensitive to the small substructures in small clusters 
than is the proportional Rap- In the constant Rap case, the mean values and the standard 
deviations of the power ratios agree with Valdarnini et al. (1999). 

From 2; = 0.5 to 2; = the power ratios become smaller in both models. This is due to 
the fact that the fraction of relaxed clusters increases with time. Like the other indicators, 
the power ratios in ACDM are larger than those in OCDM at the same redshift as expected. 

Table 9 shows the results of the KS-tests for the power ratios for the column density, 
and Table 10 shows those for the X-ray surface brightness. The flgures with asterisks in 
these tables mean that they are under the signiflcance level (5%). Figures 4-6 show the 
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Model 


Rap 


log(P2/Po) 


logiPs/Po) 


\og{Pi/Po) 


ACDM 










(^ = 0) 


0.4r2oo 


-2.90(±0.66) 


-4.38(±0.69) 


-4.89(±0.75) 




0.5r2oo 


-2.94(±0.63) 


-4.52(±0.71) 


-4.91(±0.66) 




0.6r2oo 


-2.98(±0.61) 


-4.61(±0.68) 


-4.96(±0.66) 




0.8r2oo 


-3.05(±0.66) 


-4.63(±0.69) 


-5.05(±0.72) 




1.0r2oo 


-3.15(±0.68) 


-4.70(±0.67) 


-5.20(±0.74) 


OCDM 










(z^O) 


0.4r2oo 


-3.03(±0.63) 


-4.58(±0.79) 


-5.11(±0.81) 




0.5r2oo 


-3.09(±0.68) 


-4.63(±0.74) 


-5.11(±0.80) 




0.6r2oo 


-3.15(±0.76) 


-4.61(±0.76) 


-5.10(±0.84) 




0.8r2oo 


-3.20(±0.65) 


-4.72(±0.75) 


-5.22(±0.80) 




1-Or2oo 


-3.29(±0.65) 


-4.85(±0.79) 


-5.36(±0.76) 


ACDM 










{z = 0.5) 


0.4r2oo 


-2.93(±0.69) 


-4.33(±0.67) 


-4.88(±0.69) 




0.5r2oo 


-2.92(±0.51) 


-4.37(±0.58) 


-4.89(±0.65) 




0.6r2oo 


-2.92(±0.50) 


-4.40(±0.57) 


-4.85(±0.66) 




0.8r2oo 


-2.91(±0.51) 


-4.35(±0.51) 


-4.85(±0.67) 




1-Or2oo 


-2.93(±0.51) 


-4.39(±0.63) 


-4.81(±0.64) 


OCDM 










{z = 0.5) 


0.4r2oo 


-2.76(±0.52) 


-4.29(±0.93) 


-4.78(±0.86) 




0.5r2oo 


-2.78(±0.52) 


-4.38(±0.93) 


-4.75(±0.80) 




0.6r2oo 


-2.80(±0.54) 


-4.40(±0.90) 


-4.79(±0.82) 




0.8r2oo 


-2.87(±0.61) 


-4.43(±0.77) 


-4.90(±0.96) 




1.0r2oo 


-3.02(±0.64) 


-4.44(±0.73) 


-5.02(±0.94) 



Table 6: Mean values of \og{Pml Pq) for column density 



Note. The first column describes cosmological models and redshifts. The upper two blocks are the 
data at z = and the lower two blocks are at z = 0.5. The Second column is aperture radius. The third to 
fifth columns are the mean values of power ratios of second, third, and fourth order, respectively. 
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Model 


Rap 


log(P2/Po) 


log(P3/Po) 


log(P4/Po) 


ACDM 










(^ = 0) 


0.4r2oo 


-3.92(±0.96) 


-5.86(±1.14) 


-6.52(±1.31) 




0.5r2oo 


-4.19(±0.99) 


-6.16(±1.14) 


-6.81(±1.32) 




0.6r2oo 


-4.40(±1.01) 


-6.36(±1.22) 


-7.03(±1.34) 




0.8r2oo 


-4.75(±1.05) 


-6.66(±1.32) 


-7.37(±1.41) 




1-Or2oo 


-5.04(±1.09) 


-6.93(±1.40) 


-7.67(±1.49) 


OCDM 










(z^O) 


0.4r2oo 


-4.40(±0.97) 


-6.45(±1.29) 


-7.04(±1.42) 




0.5r2oo 


-4.68(±1.03) 


-6.70(±1.39) 


-7.37(±1.59) 




0.6r2oo 


-4.90(±1.05) 


-6.83(±1.40) 


-7.51(±1.58) 




0.8r2oo 


-5.23(±1.12) 


-7.11(±1.55) 


-7.84(±1.70) 




l-Or2oo 


-5.52(±1.15) 


-7.38(±1.60) 


-8.19(±1.76) 


ACDM 










{z = 0.5) 


0.4r2oo 


-3.54(±0.83) 


-5.41(±1.08) 


-5.99(±1.10) 




0.5r2oo 


-3.79(±0.90) 


-5.64(±1.07) 


-6.30(±1.15) 




0.6r2oo 


-3.97(±0.90) 


-5.79(±1.03) 


-6.50(±1.21) 




0.8r2oo 


-4.28(±0.97) 


-6.04(±1.05) 


-6.74(±1.28) 




1-Or2oo 


-4.50(±1.11) 


-6.08(±1.31) 


-6.76(±1.49) 


OCDM 










{z = 0.5) 


0.4r2oo 


-3.63(±1.16) 


-5.35(±1.53) 


-6.03(±1.64) 




0.5r2oo 


-3.87(±1.20) 


-5.59(±1.59) 


-6.26(±1.75) 




0.6r2oo 


-4.04(±1.21) 


-5.83(±1.66) 


-6.45(±1.79) 




0.8r2oo 


-4.29(±1.35) 


-6.05(±1.77) 


-6.71(±2.05) 




1-Or2oo 


-4.59(±1.40) 


-6.30(±1.87) 


-6.98(±2.08) 



Table 7: Mean values of log(P,„/Po) for X-ray surface brightness 
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Model 




Rap 


log(P2/Po) 


log(P3/Po) 


log(P4/Po) 


ACDM 












(^ = 0) 


0.4/1- 


-^Mpc 


-5.65(±0.85) 


-7.66(±1.14) 


-8.22(±1.23) 




0.6/1- 


-^Mpc 


-6.03(±0.93) 


-7.97(±1.23) 


-8.62(±1.32) 




0.8/1- 


"^Mpc 


-6.39(±0.98) 


-8.39(±1.14) 


-8.98(±1.31) 




l.Oh- 


-^Mpc 


-6.68(±1.02) 


-8.63(±1.24) 


-9.26(±1.34) 




1.2/1- 


-^Mpc 


-6.90(±1.05) 


-8.81(±1.31) 


-9.52(±1.45) 


OCDM 












(z^O) 


0.4/1- 


-^Mpc 


-5.96(±0.89) 


-8.05(±1.20) 


-8.57(±1.29) 




0.6/1- 


-^Mpc 


-6.47(±0.96) 


-8.55(±1.21) 


-9.09(±1.40) 




o.s/i- 


-^Mpc 


-6.83(±1.05) 


-8.85(±1.39) 


-9.48(±1.58) 




i.o/i- 


-^Mpc 


-7.12(±1.08) 


-9.07(±1.47) 


-9.74(±1.64) 




1.2/i- 


-^Mpc 


-7.35(±1.12) 


-9.25(±1.56) 


-9.95(±1.70) 


ACDM 












(z = 0.5) 


0.4/1- 


"^Mpc 


-5.35(±1.04) 


-7.11(±1.46) 


-7.70(±1.43) 




0.6/1- 


-^Mpc 


-5.78(±1.20) 


-7.57(±1.54) 


-8.22(±1.69) 




0.8/?/ 


-^Mpc 


-6.13(±1.28) 


-7.99(±1.40) 


-8.58(±1.80) 




1.0/?- 


-^Mpc 


-6.36(±1.33) 


-8.17(±1.48) 


-8.68(±1.81) 




1.2/1- 


-^Mpc 


-6.58(±1.43) 


-8.15(±1.62) 


-8.80(±1.91) 


OCDM 












{z = 0.5) 


0.4/1- 


"^Mpc 


-5.58(±1.21) 


-7.28(±1.64) 


-7.95(±1.74) 




0.6/1- 


"^Mpc 


-6.05(±1.26) 


-7.89(±1.59) 


-8.48(±1.85) 




0.8/1- 


-^Mpc 


-6.39(±1.35) 


-8.12(±1.75) 


-8.74(±2.03) 




i.o/i- 


-^Mpc 


-6.70(±1.38) 


-8.44(±1.78) 


-9.14(±2.12) 




1.2/1- 


-^Mpc 


-6.94(±1.39) 


-8.67(±1.80) 


-9.34(±1.99) 



Table 8: Mean values of log(P^/Po) with constant Rap for X-ray surface brightness. 
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results of the KS-tests summarized in Tables 9 and 10. At 2; = 0, we can distinguish the two 
cosmological models by the power ratios for Ex for all Rap- Using the power ratios for a for 
most Rap except for P3/P0, we can also distinguish the two cosmological models. 

4. DISCUSSION AND CONCLUSIONS 

We have investigated morphologies of the galaxy clusters in ACDM and OCDM at 
z — and z — 0.5 using large hydrodynamic simulations. For clusters in each model we 
have calculated the axial ratios, Mint, center shifts, and multipole moment power ratios as 
statistical indicators that quantify the irregularity of cluster morphologies. 

The power ratios and the center shifts are calculated for the projected density (a) as 
well as the X-ray surface brightness (Ex). For cr, both indicators show a larger value than 
those for Ex, because the relaxation time scale of the collisionless particles is much longer 
than that of the coUisional gas particles (Frenk et al. 1999; Valdarnini et al. 1999). 

At 2; = all mean values of statistical indicators in ACDM show larger values than 
those in OCDM. These large values, which indicate large irregularity of the clusters, suggest 
more recent formation of the clusters in ACDM as expected from the analytic prediction 
(Richstone et al. 1992) and the previous numerical studies (Mohr et al. 1995; Crone et al. 
1996; Buote & Xu 1997). 

We use KS-tests to estimate the ability of the indicators to distinguish between two 
cosmological models. From the results of these KS-tests, the distributions of axial ratios in 
the two cosmological models are indistinguishable. The distributions of Mint in ACDM and 
OCDM are also similar. Using the center shifts and the power ratios for Ex at z = 0, we 
can distinguish the two cosmological models. It is possible to discriminate between ACDM 
and OCDM using the center shifts, P2/P0: and P4/P0 for a ai z — 0. At 2; = 0.5 we can 
distinguish the two cosmological models by the power ratios for Ex, but not by the center 
shifts for Ex or both the power ratios and the center shifts for a. 

Using the power ratios for Ex we can distinguish between ACDM and OCDM better 
than using those for a as shown in Table 9 and 10. This may be due to the fact that the 
relaxation time scale of dark matter that dominates a is longer than that of the gas on which 
Ex mainly depends, as described above. Since a does not settle rapidly after the cluster 
formation epoch, the power ratios for a are relatively insensitive to the difference of the 
cluster formation epochs. 

We also note that the power ratios for Ex depend more on irregularities in high density 
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Pm/ Pq 




Pap 


^ = 


^ = 0.5 


P2/P0 












0, 


■4r2oo 


2.43e-02* 


4.45e-01 




0, 


■5r2oo 


9.03e-02 


4.57e-01 




0, 


■6r2oo 


3.30^02* 


3.21^01 




0, 


■8r2oo 


1.38e-02* 


6.84e-01 




1, 


■0r2oo 


5.79e-03* 


1.04e-01 


P3/P0 












0, 


•4r2oo 


4.38C-02* 


1.16e-01 




0, 


■5r2oo 


1.23e-01 


2.20e-01 




0, 


■6r2oo 


3.37e-01 


8.40e-02 




0, 


■8r2oo 


4.28e-02* 


4.87e-02* 




1, 


■0r2oo 


1.24e-01 


2.32e-01 


P4/P0 












0, 


■4r2oo 


3.25C-02* 


4.25e-01 




0, 


•5r2oo 


1.40C-03* 


3.60e-01 




0, 


•6r2oo 


1.53e-02* 


2.81e-01 




0, 


•8r2oo 


1.43e-02* 


1.36e-01 




1, 


■Or 200 


8.48e-02 


8.97e-02 



Table 9: The results of the KS-tests for power ratios for the column density. 
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Pm/ Pq 


Pap 


^ = 


z = 0.5 


P2/P0 










0.4r2oo 


4.19e-08* 


1.03e-02* 




0.5r2oo 


7.20e-08* 


1.47e-02* 




0.6r2oo 


1.33e-06* 


2.70e-02* 




0.8r2oo 


1.47C-06* 


2.99e-02* 




1.0r2oo 


2.56e-06* 


3.30e-02* 


P3/P0 










0.4r2oo 


6.75e-06* 


2.12e-01 




0.5r2oo 


3.51e-05* 


1.42e-01 




0.6r2oo 


4.43e-05* 


9.17e-03* 




0.8r2oo 


2.24e-03* 


1.15e-03* 




1-Or2oo 


7.99e-04* 


2.43e-03* 


Pa/Po 










0.4r2oo 


2.06e-06* 


1.79e-01 




0.5r2oo 


6.83e-05* 


9.77e-02 




0.6r2oo 


4.72C-05* 


3.22C-02* 




0.8r2oo 


2.52e-04* 


1.28e-02* 




1- Or 200 


7.56e-05* 


2.68e-03* 



Table 10: The results of the KS-tests for power ratios for the X-ray surface brightness. 
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Fig. 4.— The same as Fig.3 but for PsZ-Po- 



-25- 



CD 



CO 



CO 

o 

CO 

o 

o 

CV2 
O 



"1 ^ \ 

(j{z — 0) * 
ct(z = 0.5) o 




0.4 



0.6 



0.8 



ap 



CO 
o 

-I-' ZO 

in 

o 
o 



CV2 
O 




0.4 



0.6 



2x(z = 0) * 
Sj^(z = 0.5) o. 



0.8 



ap 



Fig. 5.— The same as Fig.3 but for Pg/Po- 
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Fig. 6.— The same as Fig.3 but for P^IPq. 
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regions (i.e. in the central region) than those for a, because Ex is proportional to square of 
gas density (and square root of temperature) as p'^^T^^'^. Since small substructures which 
appear in the central region are expected to be erased more easily than large scale substruc- 
tures, the power ratios for Sx are probably more sensitive to the cluster formation epoch 
than those for a. In order to clarify this point we calculate power ratios for a"^, and then 
perform KS-tests. The results for cr^ are better than those for a, but worse than those for Ex 
reflecting the difference in relaxation time-scale. Prom these facts, we suggest that indicators 
that focus on the central regions of clusters are suitable to distinguish between ACDM and 
OCDM. 

Some previous studies (Jing ct al. 1995; Crone et al. 1996; Buote & Xu 1997) compared 
ACDM and OCDM by using pure N-body simulations in which they assume that Pgas is 
proportional to pdm- For clusters which are not well relaxed, however, this assumption may 
not be vahd. Using the power ratios for Ex Buote & Xu (1997) could not distinguish ACDM 
and OCDM, while we can. On the other hand, Jing et al. (1995) and Crone et al. (1996) 
showed that ACDM and OCDM are distinguishable by using center shifts, and we also reach 
the same conclusion. These results suggest that the power ratios are more sensitive to the 
existence of gas than the center shifts. 

We shall discuss differences in the properties of these two indicators. The center shifts 
show a large value when there are substructures whose sizes are large enough to shift a center 
of contour from the inner contour to the outer one. Since the center shifts are sensitive to 
large size substructures that can be detected in the distributions of both gas and dark 
matter in the cluster, the center shifts for the simulations with and without gas are expected 
to be similar. In contrast, the power ratios are more sensitive to small size substructures 
than the center shifts and such small size substructures are sometimes detected in only the 
dark matter distribution since corresponding gas structures are already relaxed. Hence the 
difference between the power ratios for the simulations with and without gas is expected to 
be larger than that of the center shifts. In order to confirm this explanation, we calculate 
correlation coefficients of the power ratios for Ex and a, and those of center shifts for Ex and 
(7. As expected, the correlation coefficients of the power ratios tend to be smaller than those 
of the center shifts. Therefore it is important to include a gas component in simulations in 
order to distinguish the two cosmological models using the power ratios. 

Although some authors (Evrard et al. 1993; Mohr et al. 1995) employed gas particles, 
the number of clusters was too small (8 clusters for each model). We use much larger size 
simulations (Lbox = 150/i~^Mpc) than the previous studies in order to obtain a sufficient 
number of clusters to perform statistical tests (70-80 clusters at z=0). The KS-test is 
significant when the number of samples is larger than 20 (Press et al. 1988). 
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Valdarnini et al. (1999) carried out hydrodynamic simulations for a sufficient number 
of clusters to perform statistical tests. By comparing their simulated clusters with ROSAT 
clusters, they showed that power ratios for gas in SCDM universes are in much better agree- 
ment with the ROSAT data than those for ACDM. Although this result is quite interesting, 
in order to give more conclusive results we should wait Chandra and XMM data which will 
be high resolution and high S/N and much higher resolution simulations with additional 
important physical processes such as radiative cooling, SNe feedback, and so on. 

By comparison with the power ratios in ACDM obtained by Valdarnini et al. (1999), 
those obtained in this paper arc systematically large. This difference results from the fact 
that the Rap for our power ratios is proportional to r2oo, while Rap in their paper was constant. 
Thus, our power ratios are more sensitive to small-scale irregularity for small clusters than 
theirs. Mean values of power ratios and its standard deviations in our simulations agree well 
with their results for the constant Rap case. 

Schuecker et al. (2001) showed that the fraction of clusters of galaxies with substructures 
is almost 50% for REFLEX-I-BCS clusters after correction for several systematic effects. They 
also suggested that the fraction of clusters with substructures depends on the number density 
of the clusters. Then this large fraction may be explained by the fact that a large fraction 
(~ 54% for z < 0.1) of rich clusters appear to belong to supcrclustcrs (Bahcall & Soneira 
1984). On the other hand, the size of simulation box of our and previous simulations may not 
be large enough for the realization of large supcrclustcrs, such as the Shaplcy supercluster, 
since the Shaplcy supercluster is expected to be formed from a 3.5a perturbation with a 
size of ~30Mpc (Ettori, Fabian, & White 1997). This may be the reason why the previous 
simulations of open and flat A universes do not agree with the ROSAT X-ray observation 
data. We will simulate initial conditions corresponding to the local universe to examine this 
possibility. 

Our conclusion is as follows. The power ratios and the center shifts for Sx can be 
used to distinguish between ACDM and OCDM. Since the formation histories of clusters 
in ACDM are similar to those in OCDM compared with the high-density universe, these 
indicators will be powerful probes to test the cosmological models in the non-linear regime. 

We are grateful to Vincent Eke for helpful advice. We also thank the anonymous referee 
for valuable comments. Numerical computation in this work was carried out on the SGI 
Origin 2000 at the Division of Physics, Graduate School of Science, Hokkaido University, and 
on the VPP300/16R and VPP5000 at the Astronomical Data Analysis Center of National 
Astronomical Observatory, Japan. 
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